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ON  THE  EXISTENCE  OF  SIMPLE  LIAPUNOV 
FUNCTIONS  FOR  LINEAR  RETARDED  DIFFERENCE 
DIFFERENTIAL  EQUATIONS 

by 

L.A.V.  Carvalho,  E.F.  Infante,  J.  A.  Walker 


Liapunov  functions  of  simple  form  have  been  used  for 
the  study  of  stability  properties  of  difference-differential 
equations.  In  this  paper  we  provide  necessary  and  sufficient 
conditions  for  the  existence  of  such  functions. 
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In  this  paper  we  consider,  for  the  retarded  linear  autono- 
mous difference-differential  equation 

m 

y(t)  = Ay (t)  + l Bky(t-rk),  t > 0, 
k — 1 

the  question  of  existence  of  Liapunov  functions  of  the  two  forms 

m r 0 

V (x)  = y*Ry  + l V*(0)QvV(0)d0, 

k = l ' “ T,  K 


W(x)  = sup  V* ( 0) RV (0)  . 

maxrk_<0£O 


It  is  noted  that  the  existence  of  a Liapunov  function  of 
either  of  the  foregoing  forms  is  independent  of  the  delays 
hence,  so  will  be  any  stability  conclusions  drawn  from  the  use  of 
such  functions.  This  remark  suggests  that  the  class  of  functional 
equations  of  the  above  form  for  which  such  functions  exist  is 
narrow.  In  the  paper  we  give  a characterization  of  this  class; 
moreover,  we  show  by  counterexample  that  there  do  exist  equations 
of  this  form  that  yield  asymptotic  stability  for  all  delays  {tkJksl 


and  yet  there  exist  no  Liapunov  function  of  the  above  forms  that 
can  be  used  to  establish  this  fact. 


r- 


\ 


Finally,  we  use  Liapunov  functions  of  the  above  type  to 
provide  estimates  of  exponential  decay. 
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1 n t roduc  t i on 


l . 

In  till*'  pap  or  we  cons  ulor  t ho  existence  oi  two  simple 
tvpos  of  Liapunov  functions  for  retarded  linear  autonomous 
difference-differential  equations.  The  study  of  the  stability 
properties  of  those  equations  was  pioneered  by  Yoshizawa  [1,^], 
Kazumikhin  13],  Krasovski  i |4|,  Driver  15],  and  Halo  It*]. 

Restriction  to  liapunov  functions  of  simple  form  loads 
to  stability  criteria  which  are  ovorlv  sufficient;  on  the  other 
hand,  such  functions  are  particularlv  desirable  in  specific 
applications  where  ease  of  computation  is  of  primary  importance. 
Here  we  provide  necessarv  and  sufficient  conditions  for  the 
existence  of  useful  Liapunov  functions  having  standard  simple 
forms  for  certain  difference-differential  equations. 

consider  the  di f ference  di f ferent ial  equation 

. ft 

ytt)  ■ Av(t)  ♦ ^ H.vft-i.l,  t ■*  0,  1 1 . 11 

k-1  k K 

with  initial  data 

y 1 0)  ■ x"  y ( 0) , ‘ T 1 ® 1 0 , 

where  1 ^ 1 11  • 1 max  i^,  and  the  n w n complex  matrices 

k 

A.R,  are  given,  as  is  the  initial  function  v ^ : l - 1 , 0 ] » Cn. 

Let  us  replace  the  above  equation  bv  an  abstract  evolution 


equation  on  a Banach  space  x» 


x(t)  - Fx(t)  , t > 0, 

x (0)  =x0  = (vQ  (0)  , vQ)  €£/(F)  C X.  (1.2) 


where  x = (x=(v,v)  G Cn  * V>'(  l - x , 0 ) ; Cn)  | >"*v(0) } with 


I I x | | = sup  | I v ( 0 ) | | ; the  operator  F lias  domain 

X -t^0<0  C 

,,  m 

V(I:)  = (x=  (>'  ,v)  6 X I v ' G t.(|-t,0];f  ),v'(0)«Ay  + £ Bj.v(-Tj.)},  whe 


v'(9)  = v(0),  and  is  defined  by  Fx  = (Ay  + j Rj.v(-Tj.),v') 

k 1 

for  all  x = (y,v)  G £/(F) . It  is  well  known  [7]  that  F is  the 
infinitesimal  generator  of  a C0*semigroup  {S(t))t>g  of  bounded 
linear  operators  S(t):x  X!  in  particular,  for  every 
Xq  G £/(F)  , S(-)Xg  is  the  unique  strong  solution  of  (1.2). 

We  are  interested  in  obtaining  stability  results  for  (1.2) 
(hence  (1.1))  through  the  use  of  Liapunov  functions. 


Definition  1.1:  A continuous  function  V:x  -*  fl\  is  said  to  be  a 
continuous  Liapunov  function  for  (S(t)  )t>o  on  X if  V(x)  <■  0 
for  all  x £ X»  where  V : x -*■  tR  is  given  by 


V(x)  = lint  inf  i[V(S(t)x) -V(x)  ] , x i 
t X 0 
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It  is  well  known  that  if  V is  a Liapunov  function  on  x»  then 
V(S(t)x)  <_  V(x)  for  all  t ^ 0,  x £ x • The  major  difficulty 
encountered  in  using  such  a conclusion  to  derive  stability  results 
is  the  construction  of  a suitable  Liapunov  function. 

In  [8,9]  very  complicated  Liapunov  functions  have  been 
shown  to  yield  necessary  and  sufficient  conditions  for  stability 
and  asymptotic  stability  for  (1.2), hence  for  (1.1).  On  the  other 
hand,  most  attention  [4,6,7]  has  been  centered  on  the  use  of 
functions  having  very  simple  structure;  in  particular,  functions 
of  the  form 

m f 0 

Vfx)  = y*Ry  * I v*(6)Qkv(6)d8,  x = (y,v)  e x.  Cl. 3) 

k-1  J-rk 

or  of  the  form 

W(x)  = sup  v*(0)Rv(0) , x = (y,v)  £ X.  (1.4) 

-t<0<O 

where  y*  denotes  the  conjugate  transpose  of  the  column  matrix 
y,  and  R,Qj,,  are  n * n Hermitian  matrices.  Computing  v:X  + IR 
and  iv:  X R according  to  (1.2),  we  obtain 

. m m 

V(x)  = y * (RA+ A*R+  l + 2Re  I y*RBkv^'Tk^ 

k=l  k=l 

m 

• kI1V*C-T|()Qkv(-T|(), 


Cl. 5) 
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and 


W(x)  < 


0 if  y*ny  < W(x), 


m 

maxtO,2Rey*R(Ay+  £ B.v(-t.))} 
k=l  K K 


if  y *Ry 


(1.6) 

W(x), 


for  all  x * (y,v)  G X;  hence,  V is  a Liapunov  function  on  \ 
if  and  only  if 


m 

RA*A*R*  l Q.  RB,  RB, 

k=l  k 1 


B*R  -Qj  0 . . . 0 


(1.7) 


* 

B R 
m 


0 


(i.e.,  this  matrix  is  negative-semidefinite) , while  W is  a 
Liapunov  function  on  x»  if  and  only  if 

* m 

Re:0  R(A20  + V Bk  = k)  < 0 (1.8) 

K=  1 

for  all  Iq.Zj e Cn  such  that  z*  Rz0  > z£  Rzk  for  all 

k = l,2,...,m.  We  remark  that  under  condition  (l.S)  the  Liapunov 


It  is  noted  that  he  existence  of  a Liapunov  function  of 


either  of  the  foregoing  forms  is  independent  of  the  delays 


hence  so  will  be  anv  stability  conclusions  drawn  from 


the  use  of  such  functions.  This  remark  suggests  that  the  class 


of  equations  of  the  form  (1.1)  for  which  such  functions  exist  is 
narrow.  In  the  sequel  we  give  a characterization  of  this  class; 


moreover,  we  show  by  counterexample  that  there  do  exist  equations 
of  the  form  (1.1)  that  yield  asymptotic  stability  for  all  delays 
(Tk}J=1  and  yet  there  exists  no  simple  Liapunov  function,  having 
either  of  the  foregoing  forms,  that  can  be  used  to  establish  this 
f act . 


1 1 . The  Single  Delay  Case 

In  this  section  we  consider  the  case  where  equation  (1.2) 
involves  only  a single  delay  (m  = 1 . x = i^). 

Essential  to  our  arguments  is  the  following  result. 

Lemma  2.1:  Given  n * n complex  matrices  M , N , P , wi th  M and 

N Hermitian,  the  following  two  statements  are  equivalent: 

(i)  M + N + e^P  + e'1<}>P*  < 0 fcrr  all_  E R. 

(ii)  There  exists  a Hermitian  matrix  Q such  that 


Proof : Noting  that  (ii)  implies 


[ I , e 


i4>I] 


M + Q P 

I 

P*  N - Q_ 

eiq>L 

M + N + e^P  + e‘1<J>P*  < 0 


for  all  <(>  E IR,  we  see  that  (ii)  implies  (i)  . 

A 

To  show  the  converse,  note  from  (i)  that  if  u(i<J>)  denotes 
the  fourier  transform  of  an  arbitrary  u(t),  u E (IR; Cn)  , then 


; u*(i4>)  [M+N+e 

J - 00 

equation , 


1(^P+e 


l4*P* ] u (id>) d<{>  < 0; 


hence,  by  Parseval's 


19 
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[u*(t)  (M+N)u(t)  + 2IReu*  (t)  Pu  (t- 1)  ] dt  < 0. 


Let  u(t)  = 0 for  t < 0 and  u(t)  = uk  for  k <_  t < k + 1 , 
k = 0,1,2,...;  it  then  follows  that 


00  ^ ^ 

1 [uk(M+N)uk  + 2IReukPuk_1]  £ 0 
k - 0 


for  all  uk E Cn  such  that  {uk>k=0  6 * u_i  = °- 

Define  the  functional 


* * 

J(u0)  = SUD  l [uk(M+N)uk  + 2(ReukPuk_1]  , uQ  E Cn. 

r 1°°  k=l 

{uk>k-i 


We  note  that  -Uq(M+N)Uq  _>  J(Uq)  _>  0 and,  therefore,  J(0)  = 0 

and  J(Uq)  is  continuous  at  Uq  = 0 ; furthermore,  it  is  easily 

verified  that  J(yug)  = |y|  J(uQ)  for  all  y E C. 

We  first  wish  to  show  that  is  a seminorm  on  In; 

1/2 

it  only  remains  to  be  shown  that  J (•)  satisfies  the  triangle 

inequality.  To  this  end,  consider  positive  real  numbers  a,B, 

2 2 n 

such  that  a + B =1,  and  notice  that,  for  arbitrary  w,z  E In, 


J(w  + z)  = sup  (u.  (M+N)u,  + 2IReu.  P(w+z)  + l [u.  (N+M)u, +2lReu.  Pu,  .]} 

11  k-2  K K k k-1 

Uk  k=l 


sup  {(aUj)*(M+N)(au^)+  2(Re  (uu^)  *l,(a’  w) 


l [ (“u^)  * (M+N)  (auj. ) + 2Ke  (au^)  *P  (uu^_  j ) ] } + 
ks  2 


sup  {(Bu1)*(M^N)(Bu1)  ♦ 2Re(Bu1)*P(B‘1z)  ♦ 


{uk}k=l 


l l(8uk)*(M*N)  (6uk)  ♦ 2IRe(euk)*P(3uk_1)]} 
1c  — 2 


= J(a_1w)  + J(6_1z)  = a"2J(w)  ♦ g‘2J(z). 


If  J(w)  = 0 = J ( z ) , it  follows  that  J(w  + z)  = 0;  if  not,  let 
a2  = ^2  (w)  / [ J ^2  (w)  + J1/2(z)]  , 82  = 1 - a",  arul  note  that  we 

obtain  J(w+z)  < J(w)  + 2J1^2  (w)  J1^2  (z)  + J(z).  Therefore,  J1/,2(.) 


is  a seminorm  on  (L' 


We  now  wish  to  show  that  this  seminorm  satisfies  the 
parallelogram  law.  Let  w0,zQ£  Cn , e > 0,  and  note  that  tliere 
exist  sequences  ^wk^k-l*  ^zk^k~l  ^ ^2  such  that 


9 


J(w0)  1 E + l [wk(M+N)wk  + 2IRewkPwk  l]  , 
k “ 1 

JU„)  < c + J [zk(M+N)zk  ♦ 2ReikPzk.1]; 
consequently , 

oo 

J(w0)  + J (zq)  i 2e  + j l [ (wk  + zk) *(M+N) (wk  + zk)  + 

k~  1 

+ 2«e(wk+zk)*P(wk.1+zk.1)] 

CO 

+ 7 5;  t(wk-zk)*(M+N)(wk-zk)  + 

k ~ 1 

* 2Re(wk-zk)*P(«k.1-zk.1)] 

i 2e  * Z^Ci-o^o)  * iJ("o-:o)- 

As  E > 0 is  arbitrary,  2 J ( wQ)  + 2J(zQ)  < J(wQ+z0)  + J(wp-z0); 
by  an  obvious  change  of  variables  and  quadratic  homogeneity,  the 
reverse  inequality  holds  as  well.  Hence  satisfies  the 

parallelogram  law. 

Define  q : In  x Cn  -*■  t by 

q(w,z)  = j[J(w+z)  - J(w-z)]  + ^-[J (w+iz)  - J(w-iz)],  w,z  e Cn. 

We  note  that  q is  continuous  with  q(w,w)  = J(w)  _>  0,  q(w,z)  = 
q ( z ,w) , and  a standard  argument  [10]  shows  that 
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q (ow+w  , z)  . q^,,)  foraU  „ g £ 

Consequently,  Kiesz'  representation  theorem  implies  the  existence 
of  a unique  Hcrmitian  matrix  « such  that  z.Qw  . q(w>2)  . Z.N„ 
Z,w  6 Cn. 


inallv , we  note  that,  for  every  Uj  £ £n 

qfUl,v ' J(Ui)  * {u  * ^v>Vl, 

tUk'k=2 

and  therefore,  for  every  u^Uj  <=  Cn, 


‘KOj.u,)  . Uj  (M+N)Uj  r HHeUji’Uy  < J(U(),  . q(u0,u0); 


consequently , 


U‘(Q'N)“1  * VM*N)u,  * 2ReUjPu0  < u0*tQ.N)u0, 
and  we  conclude  that  (i)  implies  (ii). 

■ 

ou^h  the  use  of  this  lemma  we  can  provide  existence 
criteria  for  Liapunov  functions  of  the  form  (1.3)  or  (1.4)  when 
equation  (1.1)  involves  only  one  delay;  i.e.. 


X(t)  = Ay  ( t)  + By(t--r),  t > 0 


(2.1) 


-a 


11 


with  m = 1,  = t,  Bj  = B.  In  this  case  the  function  (1.3)  has 

the  form 

° 

V(x)  = y*Ry  + v*  (0)Qv(0)d0 , x=(y,v)  G x»  (2.2) 

- T 

for  some  n x n Hermitian  matrices  R and  Q.  The  following 
theorem  gives  necessary  and  sufficient  condition  for  (2.2)  to  be 
a Liapunov  function  on  x f°r  (2.1). 


— 
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Proojf:  By  direct  computation  of  V (see  1.7),  we  find  that 
V(x)  satisfies  the  estimate  of  (i)  if  and  only  if  the  2n  * 2n 
Hermitian  matrix 


RA+A*R+Q+C 

BR 

R*B 

D-Q_ 

is  negative  semidef ini te ; hence,  Lemma  2.1  shows  that  (i)  is 
equivalent  to  (ii) . 

lhis  theorem  leads  to  the  following  criteria  for  stability, 

asymptotic  stability  and  instability  of  the  trivial  solution 
(hence,  by  linearity,  of  all  solutions)  of  equation  (2.1). 

I heorem  2.2:  Le  t the  re  exist  n * n He  rmi t i an  matrices  C,R, 
such  that 

RA  + A*R  + C + e 1 ^RB  + e'1,J>B*R  < 0 for  all  $ € R.  (2.3) 


Then,  for  all  t > 0, 

(i)  I_f  C > 0 and  R > 0,  the  trivial  solution 

of  (2.1)  i_s  stable . 

(ii)  if  C > 0 and  R > 0,  the  trivial  solution 


of  (2.1)  is  asymptotically  stable ; 


1 3 


n 

(iii)  It'  C > 0 and  >•  R>'q  < 0 for  so  mo  Vq  fc  C , 
the  trivial  solution  of  (.2.1)  is  unstable. 


Proof:  By  Lemma  2.1,  condition  (2.3)  implies  the  existence  of 

Q * Q*  such  that 


RA+A*R*Q+C 

B*R 


RB 


-Q 


< 0; 


hence,  Q _>  0 , V(x)  j>  v*Rv,  and  V(x)  _<  -y*Cv  for  the  function 

r° 

V(x)  = y*Ry  + v*  (<tOQv($)d$,  x = ly,v)  c X- 

' ~ T 

Using  the  semigroup  notation  of  Section  1,  we  recall 

that  C _>  0 implies  V(S(t)x)  _<  V(x),  t _>  0,  x r \;  hence, 

V(x0)  2.  V(5(t)x0)  _>  y*(t)Ry(t),  t _>  0.  Consequently,  if  R N 0 

and  I I x0 I I x < 6,  there  exists  c^  * such  that 

| |y(t)  | | j:  < cq32,  t >^0;  this  implies  that  | | S(t)  x0  | | ~ _<  cQ<5  “ 

for  all  t ^ 0.  It  follows  that  (i)  has  been  proved;  moreover, 

for  R > 0,  every  motion  S(*)xl1:IR+  -*■  x has  bounded  positive 

orbit  y(xn)  = (J  s(t)x.,.  It  is  well  known  that, for  (1.2), 
u t>0 

bounded  positive  orbits  are  precompact  [7];  hence,  making  the 
stronger  assumption  that  C > 0 and  noting  that  the  largest 

positive  invariant  set  ( 1 1 J in  (x  tr  \|\'(x)  * 0)  is  (0),  the 

Invariance  Principle  [11]  yields  the  conclusion  that  S(.t)x^  -*  0 
as  t *,  for  every  x^  ^ X»  which  proves  (ii). 
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Now  considering  (iii),  with  C > 0 and  YqRYq  < 0 for 
some  Yq  £ l’n , the  Invariance  Principle  shows  that  y(x0) 
is  not  precompact  (hence,  not  bounded)  for  some  Xq  in  each 
neighborhood  of  0 £ xi  hence,  x = 0 is  unstable  and  the  proof 
is  complete.  ^ 

We  remark  that  if  only  the  hypotheses  of  (i)  are  satisfied, 
rather  than  the  stronger  hypotheses  of  (ii),  it  may  still  be 
possible  to  employ  the  Invariance  Principle  [11]  to  conclude 
asymptotic  stability  rather  than  mere  stability;  one  need  only  be 
able  to  show  that  the  largest  invariant  set  in  {x  £ x|V(x)  = 0} 
is  {0} . 

Let  us  now  consider  the  function  W of  (1.4), 

W(x)  = sup  v*(9)Rv(0) , x = ( y , v ) £ x-  (2.4) 

- 1£0_<  0 

Recall  that,  for  this  function  to  be  a Liapunov  function  for 
equation  (2.1),  it  is  required  that 

Re  ZqR(Azq  + Bz^)  <_  0 (2.5) 

n ii 

for  all  Zq.Zj  G I such  that  zoRzO  — ziRzi*  By  making  the 
particular  choice  z^  = e^z^  we  note  that  a consequence  of 

A 


14a 


condition  12.5)  is  that  condition  (2.3)  must  hold  for  C = 0. 
Hence  existence  of  a Liapunov  function  of  the  form  (2.4)  implies 
the  existence  of  a Liapunov  function  of  the  form  (2.2),  hut  not 
conversely  unless,  in  (2.2),  we  have  that  Q * uR  for  some 
i N 0 . It  follows  that  the  class  of  equations  of  the  form  (2.1) 
which  admit  "Ruzumi khin- type"  Liapunov  functions  is  a subset  of 
the  class  that  admits  Liapunov  functions  of  the  form  (2.2). 


I 
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III.  The  Multi -do  lay  Gase. 

In  this  section  we  generalize  the  results  presented  in 
the  previous  section  so  as  to  encompass  the  case  in  which  several 
delays  occur,  as  in  equation  (1.1). 

For  this  purpose,  we  first  present  a generalization  of 
Lemma  2.1. 

Lemma  3.1:  Given  n x n complex  matrices  M , N k , P k , with  M 
and  Nj.  He rm i t i an  , k = 1,2,.  . , m , the  fol  1 ow  i ng  two  statements 
are  equivalent  : 

(.i)  There  exist  liermitian  matrices  NL*  k = l,2,...m, 

m 

such  that  M < V M.  and_,  for  each  k = 1,2,...  in, 

" . “ k=l  \ * 

M + Nk  ♦ e 1 ^ Pk  + e'l,rPk  < 0 for  all  <l>  € U. 

(ii)  There  exist  Hermit  ian  mat  rices  Q*,,  k = l,2,...,m, 


such  that 

m 
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Proof : It  is  clear  that  (ii)  holds  if  and  only  if  there  exist 

Hermitian  matrices  , k = l,2,...,m,  such  that 

m m 

I H.  > H + l Q,  and 
k=l  K " k=l  K 


k = 1,2,. . . ,m. 


Defining  Mk  = »k  - Qk  and  applying  Lemma  1 for  each  k = l,2,...,m, 
it  follows  that  (i)  is  equivalent  to  (ii)  . 


As  in  the  previous  section,  we  will  use  this  Lemma  to  pro- 
vide existence  criteria  for  Liapunov  functions  of  the  form  (1.3) 
or  (1.4)  for  equation  (1.1), 

m 

y(t)  = Ay (t)  + l Bky(t-rk),  t > 0.  (3.1) 

k= 1 

We  recall  that  we  are  interested  in  a function  V hav.iog 
the  form  (1.3), 

m ,0 

V(x)  = y*Ry  + l v*(9)Qkv(0)d9,  x = (y,v)  6 X>  (3.2) 

k=1  J‘Tk 

for  some  n x n Hermitian  matrices  R and  Qk , k = l,...,m. 


Theorem  3.1:  Given  arbitrary  n * n positive -semi definite 
Hermit  ian  matrices  C , L>k  , k =*  1 , 2 , . . . , m , the  follow  inn  two 
statements  are  equivalent : 

(i)  There  exist  R = R* , Qk  * Qk , such  that  V 

given  by  (3. 21  is  a nontrivial  Liapunov  function 
for  (3.1)  on  x and 


V(x)  £ -y*Cy  - t v*(-tk)Ukv(-Tk) , x = (y,v) 


(ii)  There  exist  R = R*  ? 0 ,M  = M.  , such  that 

m 

RA  + A*R  + C < y M,  and . for  each  k = l,2,...,m, 
k-1 


Mr  ♦ Dk  + e^RBi  + e " 1 ^ R * R < 0 for  all  <J>  G IR. 


V'  - 


The  proof  of  this  theorem  follows  immediately  from  application  of 
Lemma  3.1,  using  the  same  arguments  as  in  the  proof  of  Theorem  2.1 
This  result  immediately  leads  to 


Theorem  3.2:  Let  there  exist  n x n llermitian  matrices  C .NL  , R 

m 

such  that  RA  + A*R  + C £ M.  and . for  each  k = 1,2,. . . ,m. 


Mk  ♦ eURBk  ♦ e_1^BkR  ^ 0 for  all  $ £ 


Then , for  al  l t > 0 , k 


...  J 


1, 


* • • • »in» 


(i) 

Lf 

C 

> 0 and 

R > 0, 

the  trivial 

so  1 ut  i on 

of 

(3. 

1) 

is  stable, 

(ii) 

If 

C 

> 0 and 

R > 0, 

the  trivial 

solution 

of 

(3. 

1) 

is  asymptotically 

stable . 

(iii) 

L£ 

c 

> 0 and 

>'oR>’o  < 

0 for  some 

> o e 1 » 

the  trivial  solution  of  (3.1)  is  unstab le . 


We  omit  the  proof  of  this  theorem  since  it  parallels  that 
of  Theorem  2.2.  As  for  Theorem  2.2,  we  remark  that  it  may  be 
possible  to  conclude  asymptotic  stability  in  certain  applications 
satisfying  the  hypotheses  of  (i)  but  not  those  of  (ii);  this 
merely  requires  a more  detailed  application  of  the  Invariance 
Principle  [11]. 

The  remarks  made  in  Section  II  regarding  the  existence  of 
a "Rat umikhin- type"  Liapunov  function  (i.e.,  of  the  form  (1.4)) 
apply  here  as  well;  existence  of  this  type  of  Liapunov  function 
implies  extence  of  a Liapunov  function  of  the  form  (3.2),  but  the 
converse  need  not  hold  unless  = ctj.R,  k = l,2,...,m,  for  some 

real  a^.  > 0. 

I 
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IV . A Converse  Question  on  Asymptotic  Stability. 

If  the  hypotheses  of  conclusion  (ii)  of  Theorem  3.2  are 
satisfied,  then  we  are  guaranteed  asymptotic  stability  of  all 
solutions  of 


m 

y(t)  = Ay (t)  + l B.y(t-Tk),  t > 0,  (4.1) 

k=i  k K 

irrespective  of  the  values  of  the  delays  ik,  k = l,...,m.  It 
follows  that  no  function  V of  the  form  (1.3)  will  be  a Liapunov 
function  with  R > 0,  V(x)  _<  -y*Cy,  C > 0,  unless  equation  (4.1) 
is  asymptotically  stable  for  all  delays.  Given  this  observation, 
and  recalling  the  power  of  the  Invariance  Principle  [11],  it  is 
natural  to  ask  the  following  question:  given  that  the  trivial 
solution  of  (4.1)  is  asymptotically  stable,  irrespective  of  the 
delays  ik>  k = l,...,m,  does  there  necessarily  exist  a Liapunov 
function  of  the  form  (1.3)  with  R > 0?  This  is  essentially  the 
question  posed  by  Hale  (see  [7],  p.108). 

To  investigate  this  question,  let  us  recall  a recent  result 
[12]  that  characterizes  those  (m+l)-tuples  of  real  matrices 
(A,Bj , . . . ,Bm)  such  that  the  trivial  solution  of  (4.1)  will  be 
asymptotically  stable  irrespective  of  the  delays. 

Theorem  4.1  [12]:  The  trivial  solution  of  equation  (4.1)  is. 

asymptotically  stable  for  any  set  of  xk  >.  0 if  ajid^  only  if 

m 

(i)  the  real  parts  of  the  eigenvalues  of  A + y B,, 
k = 1 K 


are 
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negative,  and 

m i <(>k 

ii)  for  all  4>,,  G |R,  A + ) e B,.  has  no  nonzero 

k=l  K 

imaginary  eigenvalue . 

Thus,  in  the  case  of  real  matrices  A,Bk,...,B  , our  question 
reduces  to  the  following:  do  the  hypotheses  of  Theorem  4.1  imply 
the  existence  of  Hermitian  matrices  R>0,  C^O,  Dk  ^ 0 and 
Mk , k = l,...,m,  such  that  statement  (ii)  in  Theorem  3.1  is  true? 

This  question  has  an  affirmative  answer  in  the  scalar 
case  (i.e.,  A = a G tR,  Bk  = bk  G (R)  . Indeed  the  conditions  of 
Theorem  4.1  become,  in  this  case 


m 


l bk  < °, 

k=  1 K 

(4.2) 

m 

I 1 bk | < 0. 

k = l 

(4.3) 

Hence,  letting  R = 1 , C = 0 , Dk  = 0 , and  = -2|bk|,  it  is 
immediately  seen  from  (4.3)  that  all  conditions  of  statement  (ii) 
in  Theorem  3.1  are  satisfied,  and,  therefore,  there  exists  a 
Liapunov  function  of  the  form 


V(x) 


m 

k=l 


Qklv(0) l2de, 
k 


x = (y,v)  G x, 


(4.4) 


for  some  . ,Qm  £ R,  such  that  V(x)  <_  0.  From  the  proof 

of  Lemma  3.1,  it  is  easily  seen  that  the  Qk  must  satisfy 


• 1 


I 


- 


L 
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2|bkl*Qk  b, 


-Qi 


£0,  k = 1,2,. . . ,m; 


hence , = |bjJ  is  the  only  possible  choice,  and  we  obtain 


• n!  7 m - 

V ( x ) = 2 (a+  l |bk  1)1x1“  - l I b,  | [y-  (sgn  b,  )v(- t.)  ] “ 
k= 1 K k= 1 K k k 


for  all  x = (y,v)  e X. 

Alt  hough  C Ji  0 and  result  (ii)  of  Theorem  3.2  may  not 
apply,  the  Invariance  Principle  [11]  can  be  employed  directly  to 
show  that,  under  condition  (4.21,  this  Liapunov  function  does 
establish  asymptotic  stability. 

The  above  question,  however,  has  a negative  answer  at  the 
vector  level,  as  we  can  show  through  a counterexample.  Indeed, 
in  (4.1)  let  in  = 1 and 


-1 

-1 

a 

cf~ 

A = 

1 

-1 

, B,  = 

0 

-a 

for  some  « 6 IR.  It  is  easily  verified  that  the  conditions  of 
Theorem  4.1  are  satisfied  if  and  only  if  «“  < 2.  On  the  other 
hand,  it  is  not  possible  to  satisfy  condition  (2.3)  with  any 
R = R*  > 0,  C = C*  > 0,  when  a1-  > 1. 


Indeed,  suppose 
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1 < a <2  and 


R(A+e1<J>B1)  + (A+e^Bj)  *R  _<  -C  for  all  4>  e IR,  (4.5) 


where  C = C*  > 0 and 


1 B 


B*  6 


Bee,  oteiR,  6 > |b| 


Satisfaction  of  (4.5)  implies  that,  for  all  <j>  e K, 


-2+2acos  4>  + 2ReB  - l-23  + 6-2aBsin 


- l-2B  + 6-2otBsin  4>  -2ReB-26-2<$acos  4> 


for  this  to  be  so,  it  is  necessary  that 


-26(l+acos  <t>)  2ReB  £ 2(l-acos  40 


for  all  4>  E IR,  and  this  is  impossible  for  1 < a4- . Hence, 

2 

for  1 ^ a <2,  we  have  an  example  displaying  asymptotic 
stability  for  all  t _>  0,  but  for  which  no  nontrivial  function 
of  the  form  (1.3)  is  a Liapunov  function. 


V.  Additional  Remarks . 

In  this  section  we  present  two  additional  results  based 
on  the  use  of  a simple  Liapunov  function  of  the  form  (1.3).  The 
first  of  these  provides  a sufficient  condition  for,  and  an 
estimate  on,  exponential  decay  of  solutions.  The  second  provides 
sufficient  conditions  for  stability  and  asymptotic  stability 
that  are  simpler  than,  but  not  as  "necessary"  as,  those  of 
Theorem  3.2. 

In  result  (ii)  of  Theorem  3.2,  and  also  (ii)  of  Theorem 
2.2,  we  have  concluded  only  asymptotic  stability;  without  addi- 
tional assumptions,  this  conclusion  can  be  strengthened  to  ex- 
ponential stability  by  employing  yet  another  function,  closely 
related  to  the  function  V of  (1.3), 

m ,0  6(G  + t,) 

V x ( x ) = y*Ry  + [ c k v*  (0)Q,  v(0)d0,  x = (y,v)  (.  x, 

k=l  '-T.  K 

(5.1) 


where  R = 
puting  v6 


R\  Qk 

x *•  'R 


* 

= , and  6 is  a positive  real 

according  to  (1.2),  we  obtain 


numbe  r . 


Com- 


m 5 t . m 

Vs(x)  - y*  ( RA+A*  R+  e KQ.)y  + 2Re  l y*RB.  v(- 1.) 

k = l k=»l 

m 

- I v‘(-Tk)QkV(-Tk) 
k=  1 

m rO  3(0+ri) 

-<S  £ e v* ( O)Q^v(O) do , x = (>’,v)  t x- 

k-ll-Tk 


(5.2) 
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Under  the  hypotheses  of  (ii)  in  Theorem  3.2  there  exist 
(.  = C*  > 0 , R = R*  > 0 , Qk  = such  that  the  Hermitian  matrix 


m 


RA+A*R+C+  l Q Rfi 


k = 1 


* 

BlR 


B2R 


•k 

BmR 


-Q 

o 


l 


l 


RB. 


-Q 


2 • 


RB 

r 

0 

0 


m 


• -q 


m 


is 


negative-semi definite.  Consequently,  Qk  > 0 and  the  use  of 


these  matrices  in  leads  to 


in  +6t, 


V6(x)  < -6  6(x)  + y*(RA+A*R+fiR+J>e  kQk)y  ♦ 2ReJ  y*RBkv(-ik) 


k=  1 


k=  1 


m 


1 v*(-Tk^kv(‘Tk) 
k=l  K K K 


m 6 T. 


i-6V6(x)  - y*[C-6R-  l (e  -l)Qk]y. 

k - 1 

Since  C > 0,  we  may  now  choose  a 6 > 0 so  small  that 
V6(x)  <.  _<$Vg(x),  which  leads  to  the  estimate 


■fit 


v<$(S(t)x)  _<  e V<.(x)  for  all  t > 0,  x e X.  (5.3) 
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Consequent ly , 


>’ * ( t ) Ry  ( t ) < e"'St  {Vq(0)Rvq(0) 


■f 


m 


l 

k=i 


<5  T , r 0 * 

c J v0(0)V(0)J0} 


for  all  t _>  0 , implying  the  existence  of  ^ 1 such  that 

I I S(t)xQ  I I ^ 1 cQe  6t  | | xQ|  | “ for  all  t _>  0,  xQ  t X«  It  follows 

that  the  hypotheses  of  (ii)  in  Theorem  3.2  (hence,  of  (ii)  in 

Theorem  2.2)  actually  imply  exponential  stability  of  the  trivial 

solution,  with  exponent - 5t/2  such  that  5 > 0 and 
m 5 t j. 

C 2 3R  * } (e  -1)Q,;  hence,  6 depends  on  the  delays  t, , 

k=l  K k 

k = 1 , 2 , . . . , m . 

Turning  to  the  second  remark  we  wish  to  make,  we  refer 
again  to  the  function  V of  (1.3),  recalling  that  V is  a 
Liapunov  function  for  (1.2)  if  and  only  if  (1.7)  holds;  i.e., 


— m 

RA+ A*R+  l Q. 


RB, 


RB, 


Q j 0 . . . 0 


< 0 


(5.5) 


0 


26 


* 

for  some  R = R* , Qk  = Qk , k = 1 ,2 , . . . ,m.  Theorem  3.1  shows 
this  condition  to  be  satisfied  if  and  only  if  statement  (ii)  of 
that  theorem  is  true  for  C = 0,  = 0,  k = 1,2 m.  Un- 

fortunately, this  statement  may  be  difficult  to  verify  and  does 
not  suggest  an  explicit  construction  for  the  matrices  Qk.  In 

certain  applications  it  may  be  desirable  instead  to  simply  "guess" 

* 

= Qk,  and  then  search  for  some  R = R* 
such  that  (5.5)  holds.  We  will  now  describe  a simple  technique 
of  this  type. 

Condition  (5.5)  implies  that  Qk  _>  0 with  null  space  no 

larger  than  that  of  RBk.  Hence,  let  us  suppose  Qk  to  be  of 
* * 

the  form  Qk  = B^RZ^RB^  for  some  > 0,  k = l,2,...,m. 

* 

Then  (5.5)  will  be  satisfied  if  there  exists  R = R*,  Zj,  = Z^  > 0, 
such  that 

m ^ 

RA  + A*R  < - l (BkRZkRBk+Yk)  (5.6) 

k_  1 

where  is  a Hermitian  matrix  that  inverts  on  the  range 

of  RBk(Yk=Zk1  if  det  | RBk | ^ 0)  . In  fact,  a careful  review  of 

the  proofs  of  Lemma  2.1  and  Lemma  3.1  reveals  that  if  (5.5)  can 

* 

be  satisfied  by  any  R = R*  f 0 , Qk  = Qk , then  it  can  be  satis- 

* 

fied  by  some  R = R*  i 0 , Qk  = Qk , with  Qk  of  the  form  we  have 

suggested.  Hence,  satisfaction  of  (5.6)  by  some  R = R*  f 0 , 

* 

Zk  = Zk  > 0 , is  necessary  and  sufficient  for  the  existence  of 
any  non-trivial  Liapunov  function  of  the  form  (1.3)  for  equation 
(1.2). 


a plausible  form  for  Qk 
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Condition  (.5.6)  certainly  will  be  satisfied  if  there 
exists  R = R*  ^ 0 and  real  numbers  6k  > 0 such  that 

m * i 

RA  ♦ A*R  ♦ l CBkRkRBk+6k  R)  < 0.  (5.7) 

k — 1 

Although  much  more  restrictive  than  our  necessary  and  sufficient 
conditions,  the  comparatively  simple  condition  (5.7)  is  certainly 

sufficient  to  ensure  that  V is  a Liapunov  function  on  x for 

* 

(1.3)  with  Qk  BkBgRBk»  and  that  the  trivial  solution  is  stable 
if  R > 0;  moreover,  V is  given  by 

. m * , 

V ( x ) = y * [ RA+A*R+  V (‘,kBkRBk  + 6k  10  ) >' 

k “ 1 

' 6k*  ly"6kBkV^"Tk^  R(y'BkBkv('Tk)  1 » 

x = (y,v)  e x. 

Hence,  under  (5.7),  if  R > 0,  it  is  possible  t'  use  the  In- 
variance Principle  111]  to  provide  conditions  sufficient  for 
asymptotic  stability;  for  example,  it  is  easily  found  that  the 
triviil  solution  is  asymptotically  stable  if  R = R*  > 0 and 
6k  > 0 , k = 1 , 2 , . . . ,m, 
strictly. 


are  such  that  inequality  (5.7)  holds 
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